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Abstract. We present the first results on the identifica- 
tion and study of very distant field galaxies in the core of 
cluster-lenses, using a selection criterium based on both 
lens modelling and photometric redshifts. We concentrate 
on two multiple-imaged sources at z = 4.05 in the clus- 
ter Abell 2390. The two objects presented in this paper, 
namely H3 (cusp arc) and H5 (fold arc), were identi- 
fied through lens modelling as multiple images of high- 
redshift sources at z>3.5 (Kneib et al. 1999). We confirm 
the excellent agreement between this identification and 
both their photometric redshifts and morphologies. Our 
CFHT/WHT program for a systematic redshift survey of 
arcs in clusters has allowed to obtain a set of spectra on 
three different images at z ~ 4: the brightest image of H3, 
whose redshift was already confirmed by Frye & Broad- 
hurst (1998), and the two brightest images of H5. The later 
is then confirmed spectroscopically as a multiple image, 
giving a strong support to the lens model. The main fea- 
ture in each of these spectra is a strong emission line, iden- 
tified as Ly a, leading to z — 4.05 for both H3 and H5. The 
spectrophotometric properties of these galaxies are stud- 
ied, in particular the degeneracy in the parameter-space 
defined by the SFR type, age, metallicity and reddening. 
H3 and H5 are intrinsically bright and clumpy galaxies 
(M B to M* B -2 magnitudes), located ~ lOO/i^ 1 kpc apart 
on the source plane, with mean metallicities compatible 
with a fraction of solar or even solar values. These results 
seem to favour a hierarchical merging scenario, where we 



Send offprint requests to: R. Pello roser@obs-mip.fr 
* Based on observations collected with the Canada-France- 
Hawaii Telescope (CFHT), the William Herschel Telescope 
(WHT) and the Hubble Space Telescope (HST). 



are seeing a relatively evolved phase in these two z ~ 4 
objects, with stars forming locally and efficiently. 

Key words: Galaxies: clustering - evolution — Cosmol- 
ogy: observations - gravitational lensing 



1. Introduction 

The identification and study of very high-redshift galax- 
ies (z > 4) is probably one of the most direct methods 
to constrain the scenarios of galaxy formation and evolu- 
tion. The observation of such galaxies presents a major 
challenge complementing the statistical studies of lower 
redshift galaxies (CFRS: Lilly et al. 1996; Hawaii Deep 
Fields: Cowie et al. 1996; HDF: Sawicki et al. 1997; Lowen- 
thal et al. 1997; Steidel et al. 1996b). The key information 
to retrieve is the star formation history in the universe 
and the evolution of the different morphological types of 
galaxies. The main problems are the identification of high- 
z galaxies and the construction of a sample as bias-free as 
possible. Steidel et al. (1992, 1993, 1996a, 1997) first used 
the Lyman dropout technique to photometrically identify 
3.0 < z < 3.5 galaxies in empty fields, and they suc- 
ceeded in confirming spectroscopically a large sample of 
them with the Keck telescope. A similar technique has 
been used on the HDF to identify 2.2<z<3.5 galaxies 
(Steidel et al. 1996b; Lowenthal et al. 1997) also con- 
firmed with the Keck telescope. The number of massive 
star-forming galaxies at high redshift can be used to con- 
strain the cosmological models (see Baugh et al. 1998). 
All these galaxies are in general very compact in their 
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rcstframc UV, with the bulk of their star formation lo- 
cated on high surface-brightness regions. Their morpholo- 
gies might give some hints on the physical processes in- 
volved in galaxy evolution. In particular, they might pro- 
vide useful constraints on the formation of the different 
systems (spheroids, bulges, disks), and its close relation 
with metal enrichment timescales (see Trager et al. 1997 
for a detailed discussion). At z ~ 3 to 4, the correspond- 
ing lookback time is 13.7 to 14.4 Gyr, 15 Gyr being the 
present age of galaxies (H = 50 km s _1 Mpc -1 , go = 0.1 
and A = 0). Thus, such high redshift systems provide a 
direct measure of the early star-formation processes. 

The photometric technique based on the identification 
of Lyman dropouts in the U-band has shown to be a useful 
tool to select 2.2 < z < 3.5 galaxies, and it has been success- 
fully extended to the B-band to locate 3.5 < z < 4.5 candi- 
dates (Steidel et al. 1998). Nevertheless, two major biases 
will appear in all such selection techniques: one towards 
intrinsically luminous galaxies because the sample is lim- 
ited in apparent magnitude, and another towards galaxies 
with strong star-formation activity, because the rest-frame 
UV will be seen at the visible wavelengths. When high- 
rcdshift galaxies are selected through photometric redshift 
techniques based on a large wavelength domain, including 
near-IR, and close to critical-lines of cluster-lenses, the re- 
sulting sample is less sensitive to these biases. This paper 
presents the first z ~ 4 results on this original method to 
build-up and study an independent sample of z > 2 galax- 
ies, by combining photometric redshifts (visible and near- 
IR) and lens modelling. One of the first examples of such 
a technique on lensed galaxies is the z = 2.51 star-forming 
object in A2218 (Ebbels et al. 1996), and a number of re- 
cent discoveries of z > 2.5 lensed-galaxies in clusters, some- 
times serendipituous, strongly encourages this approach 
(Ycc et al. 1996a; Trager et al. 1997; Franx et al. 1997; 
Soifer et al. 1998; Seitz et al. 1998; Frye & Broadhurst 
1998; Bunker et al. 1998). 

The two z ~ 4 lensed galaxies discussed in this paper 
were identified in the cluster Abell 2390 (z = 0.231), a 
gravitational lens which has been extensively studied. The 
first detailed photometric and spectroscopic surveys on 
this target were performed by Le Borgne et al. (1991), and 
recently enlarged by the CNOC group (Yee et al. 1996b; 
Abraham et al. 1996; Carlberg et al. 1996). This cluster, 
which is known as a strong X-ray emitter (Ulmer et al. 
1986; Pierre et al. 1996), has an elongated shape (e = 
(a — b)/(a + b) = 0.55, based on the luminosity map of 
the photometrically selected cluster-galaxies, Pello et al. 
1991) and a high velocity-dispersion (1090km/s, Abraham 
et al. 1996). A near-IR photometric survey is also available 
(Miralles et al. 1999, in preparation). We will focus here 
on two multiple images at z ~ 4, namely H3 and H5, 
according to the nomenclature by Kneib et al. (1999, see 
also Figure 1). The redshift of the former had already 
been confirmed by Frye & Broadhurst (1998) at the Keck 
telescope. 



A plan of the paper follows. First, we give a summary 
of the spectroscopic and photometric data in Section 2. 
Section 3 is devoted to the study of the main properties 
of H3 and H5, namely the spectroscopic and photometric 
redshifts, the lensing properties and the morphology of 
these objects. The stellar population of these objects is 
characterized in Section 4. In Section 5 these results are 
discussed, with a special attention to the implications for 
galaxy formation and evolution. Throughout this paper, 
we assume H = 50 km s _1 Mpc -1 and A = 0. 

2. Observational data 

2.1. Photometry 

The photometric data were obtained during a number of 
runs since 1988, at the 3.6m Canada-France-Hawaii Tele- 
scope, the 3.5m telescope of Calar Alto (CAHA, Spain) 
and the 2.5m Isaac Newton Telescope (Spain). We have se- 
lected the images obtained under photometric conditions. 
The total number of filter bands considered is 10 when 
we add the two F555W (V w ) and F814W (I w ) WFPC2 
images from our HST program (PI. B. Fort), although the 
two B bands have very similar central wavelengths (but 
different widths). Table 1 provides with the central wave- 
lengths and widths of the different filters. In the case of 
the 2 pairs (R and r, I and Iw), the difference between 
the central wavelengths is of the order of ~ 300A, and the 
filter widths are different. We use all these data as comple- 
mentary information in the following discussion, and we 
consider different independent sets of filters to constrain 
the properties of the stellar population (see Sect. 4). All 
magnitudes are given in the Vega system, and we use the 
SED of Vega and the detailed filter and detector responses 
to set the zero points when converting magnitudes into 
fluxes. As most of these data have been published else- 
where, we give only a summary of the main characteris- 
tics of the images in Table 1, including the detection levels 
and references. The 2a magnitude is defined here as the 
magnitude of an object with 4 connected pixels 2a above 
the sky level. 

2.2. Spectroscopy 

The spectroscopic data set used in this paper comes from 
three different runs: one at the CFHT and two at the 4.2m 
William Herschel Telescope (WHT) . H3 was one of the tar- 
gets of the spectroscopic survey performed by Bezecourt 
and Soucail (1997) at CFHT in August 1995. The spec- 
trograph used was MOS/SIS (Le Fevre et al., 1994), with 
the V150 grism, providing a low dispersion of 7.3 A/pixel. 
H5a and H5b were observed during two separate runs at 
the WHT, in September and June 1996 respectively, with 
the LDSS-2 multiobject spectrograph (Allington-Smith et 
al. 1992). The grism used was the medium blue, with a 
dispersion of 5.3 A/pixel. The total exposure times were 
6.7 ksec (2 x 2500 sec + 1700 sec) on H5b and 9.9 ksec 
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Table 1. Characteristics of the images and detection levels. 
2<r magnitudes correspond to objects with 4 connected pixels, 
each 2 a above the sky level. Observations were carried out at 
the CFHT (1), the INT (2), HST (3) and CAHA (4). 
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(2 x 3600 sec + 2700 sec) on H5a. The reduction was per- 
formed using standard IRAF procedures as well as our 
own software packages. 

3. Photometric and spectroscopic study of H3 and 
H5 

3.1. Lensing properties, spectral energy distribution and 
morphology of the sources 

Hereafter we use the new lens model of A2390 by Kneib et 
al. (1999), which is a refined version of the earlier model 
presented in Pierre et al. (1996), based on lensing and X- 
ray data, taking into account the new constraints given 
by the HST images. According to this lens modelling, H3 
and H5 are multiple images of two high-redshift sources at 
z > 3.5. This redshift estimate has been confirmed through 
photometric redshift techniques (filters from B to K) and 
spectroscopy. Figure 1 displays a zoom on the different 
components of H3 and H5, as well as their location in the 
cluster. H3a-b-c is an impressive cusp arc showing several 
bright knots which are identifiable in each different image. 
H5 is a fold arc with the two radial components showing 
similar morphologies. A third faint image is predicted for 
this source, but unfortunately it lies on the edge of the 
Planetary Camera, where the exposures have a poor S/N. 
Table 2 summarizes the photometry of all these compo- 
nents except H3c, which is too close to the bright cluster 
galaxy. The different colors have been obtained through 
aperture magnitudes computed within the same region in 
all the filters, correcting for sampling and seeing effects. 
The averaged colors of a bright E/S0 cluster galaxy are 



also given for comparison. The surface brightnesses and 
colors are compatible with the multiple-image hypothesis 
within the photometric errors (Am ~ 0.1 typically, and 
Am ~ 0.3 in the near-IR). The faint images H5a and H3b 
are surrounded by bright objects (see Fig. 1), and they are 
hardly detected in the near-IR and g where the sampling 
(> 0'.'5 per pixel) and/or the seeing conditions are poor. 
Moreover H5a and H5b are hardly detected in B. In these 
particular cases, photometric errors are at least Am ~ 0.5, 
and these magnitudes are indicated by ":" in Table 2. It 
is worth noting that the identification of multiple images 
based on the similarity of the spectral energy distributions 
(SED) is more easily obtained with extended wavelength 
coverage. 



Table 2. Photometry of the different H3 and H5 components, 
compared to the mean values for E/S0 cluster galaxies (see 
text). 
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The first spectrum of H3 obtained by Bezecourt and 
Soucail (1997) shows a single emission line which was as- 
signed to [O II]3727 A, giving z — 0.647, a value which is 
in fact incompatible with both the lens modelling and the 
photometric redshift. On the contrary, when it is correctly 
assigned to the rest-frame Lya, it gives z = 4.05. Subse- 
quent observations of H3 by Frye & Broadhurst (1998) at 
the Keck Telescope obtained a redshift of z = 4.04 through 
a higher quality spectrum. This time, the redshift is based 
on several absorption lines, the Lya centroid being slightly 
shifted redwards with respect to the absorption system, 
in good agreement with the position of the line found by 
Bezecourt and Soucail (1997). 

In the case of H5a and H5b, both spectra show a strong 
emission line (Fig. 2). When this line is identified as Lya, 
the corresponding redshifts are z = 4.049 ± 0.003 and 
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Fig. 1. Image of Abell 2390 in logarithmic scale, obtained by coadding the final HST images in Vw and Iw- A zoom on H3 
and H5 is also shown, with the identification of the different components. The critical lines at z = 1, 2.5 and 4 (from the inner 
to the outer Dart. respectively! are also displayed, according to the model bv Kneib et al. ("1999 1 ). 
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z = 4.052 ± 0.002 for the a and b components respec- 
tively. Thus, this result confirms that H5a and H5b are 
two images of the same z — 4.05 source, in perfect agree- 
ment with the lens model. When taking the Lya centroid 
to compute the redshift, the two multiple images H3 and 
H5 are at the same redshift within the errors. Hereafter we 
take the Lya-based redshift of z = 4.05 for both sources, 
because we have no estimate based on absorption lines for 
H5. The measured equivalent width of Lya is relatively 
high for H5, W\ = 273 A for the brightest image, corre- 
sponding to a rest frame value of 54 A. We obtain a rest 
frame upper limit of ~ 75 A for the faintest image, where 
the continuum is hardly detected. These values are similar 
to the ones observed by Hu et al. (1998) in their sample of 
emission-line galaxies at z ~ 3 to 6. The equivalent width 
of H3 is smaller, the upper limit obtained from the CFHT 
spectrum being ~ 33 A. 

The gravitational amplification computed for the 
brightest images H3a and H5b is 2.3 ± 0.3 magnitudes 
in both cases, and the predicted difference between the 
two images are Am a & = 0.4 and 0.8 magnitudes for H3 
and H5 respectively. The corresponding measured dif- 
ferences (Table 2) are also consistent with this model, 
Am ab (H3) = 0.4 ± 0.13 and Am ab (H5) = 0.7 ± 0.16. 
These amplification factors are surface-averaged values 
computed with the Kneib et al. (1999) lens model. 

We have applied a lens inversion procedure to restore 
the morphology of these sources on the source plane at 
z = 4.05. This method is close to the LensClean algo- 
rithm (Kochanck & Narayan 1992; Kneib et al. 1994). 
Figure 3 displays the resulting morphology and location 
of the two objects on their source plane. The true sepa- 
ration between H3 and H5 in the source plane would be 
16", corresponding to a linear separation of 150(252)/i 50 1 
kpc with qo = 0.5(0.1). Note that these sources are not 
resolved in their width, the lens inversion being limited 
by the resolution of the composite HST images (0.13" in 
Vw and Iw)- Both sources are extremely clumpy. H5 con- 
sists of an alignment of two compact and bright blobs, 
less than 0.1" apart (0.9(1.6)/i^~ kpc), and a faint ex- 
tended component of ~ 0.3" (2.8(4.7)/i^ 1 kpc), the to- 
tal length being - 0.4" (3.7(6.3)/i 50 1 kpc). H3 is more 
elongated than H5, and it displays four small and bright 
subclumps. Each one of these subclumps has less than 
0.1" of diameter, and the total length of the structure is 
about - 0.8" (7.5(12. e)^ 1 kpc). The orientation of the 
two sources is similar (see Fig. 3). Compared to the mor- 
phologies and sizes of the 2.5<z<3.5 sample by Steidel 
et al. (1996a), typically ~ 0.5" - 1" for the resolved cores 
(~ 3.6 to 7.3/igQ 1 kpc with q — 0.5; see also Giavalisco et 
al. 1996), the bright subclumps of H3 and H5 are more 
compact, but the total length of the emitting region is of 
the same order. From the morphological point of view, H3 
and H5 are similar to the 3.3 < z < 4.0 lensed sources found 
behind C10939+4713 (Trager et al. 1997). 




4500 5000 5500 6000 6500 7000 7500 

observed wavelength (A) 




900 1000 1100 1300 1300 1400 1500 
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Fig. 2. Top: Mean spectra of the b (top) and a (bottom) com- 
ponents of H5, each of them showing a strong emission line. 
Bottom: Averaged spectrum of the two H5 components, show- 
ing the emission line identified as Lya at z = 4.05. No other 
features are clearly identified on the continuum. 

3.2. The photometric redshift approach 

H3 and H5 are the first two high-redshift objects spec- 
troscopically confirmed in our sample of photometrically 
selected candidates at z > 2 in this cluster. A photo- 
metric redshift method has been used to identify such 
candidates in the cluster core, close to the critical lines. 
Photometric rcdshifts were derived according to the stan- 
dard minimization method described by Miralles, Pello & 
Roukema (1999) and Miralles (1998). The observed SED 
of each galaxy, as obtained from its multicolor photom- 
etry, is compared to a set of template spectra. The new 
Bruzual & Chariot evolutionary code (GISSEL98, Bruzual 
& Chariot, 1993, 1998) was used to build 5 different syn- 
thetic star-formation histories, each with solar metallicity 
(Zq): a burst of 0.1 Gyr, a constant star-formation rate, 
and 3 fi models (exponential-decaying SFR) with charac- 
teristic time-decays matching the present-day sequence of 
colors for E, Sa and Sc galaxies. The template database 
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Fig. 3. Restored Iw image of the two objects H3 (left) and H5 (right) on the source plane at z — 4.05, as obtained from the 
lens-inversion procedure. Both sources are clumpy, elongated and exhibit the same orientation. The distance between H3 and 
H5 is 16". At z — 4.05, 1" corresponds to a linear separation of 9. 4(15. 7)h^ kpc with go = 0.5(0.1). Isocontour plots display 
the center and profile of the different clumps, according to a linear scale which is identical for the two objects. 



includes 255 synthetic spectra. The intergalactic absorp- 
tion in the Lyman forest is modelled using the average flux 
decrements Da and Db, according to the original defini- 
tion by Oke & Korycansky (1982). Da and Db correspond 
respectively to the continuum depression between Lya and 
Ly/3, and between Ly/3 and the emission Lyman limit. The 
prescriptions for the redshift distribution of Da and Db 
are taken from Giallongo & Cristiani (1990), and they are 
in good agreement with those given by Madau (1995) in 
the common redshift domain 2.5 < z < 5. When applied to 
our data, the photometric redshift method identifies ~ 30 
sources at z > 2 in this field, most of them too faint to 
be confirmed spectroscopically using 4m telescopes. The 
high gravitational amplification of H3 and H5, and their 
strong Lya emission makes it possible in these two par- 
ticular cases. The results obtained on the whole field of 
A2390 and two other cluster-lenses will be discussed in 
more details elsewhere (see a preliminary version in Pello 
et al., 1998). 

In the case of H3 and H5, we have examined the sensi- 
tivity of our photometric redshift estimate as a function of 
the relevant parameters, namely SFR, age and metallicity 



of the stellar population. This exercise is especially impor- 
tant for H5, because the spectroscopic redshift is based on 
a single line. Figure 4 presents a likelihood map for this ob- 
ject, showing a good agreement between the spectroscopic 
(z = 4.05) and the photometric redshift (z = 3.90+Q%f, 
where the error bar corresponds to a la level). This map 
was obtained using the set of 10 different SFRs presented 
in the next section, representing templates with a range 
of metallicities. Each point on the redshift-age map corre- 
sponds to the best fit of the SED obtained across the SFR- 
metallicity space. The dark regions in Figure 4 (above 99% 
confidence level) result from the overlap of the 5 metal- 
licities considered (age vs. metallicity degeneracy). The 
redshift region around z ~ 4 appears as the most likely 
solution for this object. The map obtained when using the 
solar-metallicity set above mentioned is qualitatively the 
same. In the case of H3, we obtain a photometric redshift 
of z = 3.807o6o using the same procedure, thus in good 
agreement with the spectroscopic value, but with a larger 
error bar. 
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wavelengths seen in the visible bands. Two kinds of SFRs 
were considered: a single stellar population (instantaneous 
burst), and a continuous star- forming system, both with 
the Salpeter (1955) IMF, with upper and lower mass-cutoff 
of O.1M < m < 125Af©, and an extinction law of SMC 
type given by Prevot et al. (1984). When computing quan- 
tities related to the stellar mass involved in one of these 
regions, we take the above mass limits for star-formation, 
but only ~ 1/3 of the stellar mass corresponds to stars 
with m > IMq with this particular IMF. 



3.0 4.0 
redshift 

Fig. 4. Photometric redshift likelihood-map of H5 showing 
the excellent agreement with the spectroscopic redshift of this 
object, which is contained within the region at 68% confidence 
level. The scaling displayed at the top represents the likelihood 
value associated to the x 2 . The shaded area enclose the 3a 
contour (confidence level of 99%, or a likelihood value of 1%). 
The shaded region on the lower part of the map is excluded 
because of age-limit considerations for the stellar population 
(stars cannot be older than the age of the universe, with Ho ~ 
50 km s" 1 Mpc" 1 and q = 0). 



4. Constraining the stellar population in H3 and 
H5 

We use the SEDs of H3 and H5, determined from broad- 
band photometry, to infer the relevant parameters for the 
stellar population, in particular the total SFR. The SEDs 
of these objects can be fitted by different synthetic stellar 
populations, and there is a degeneracy to consider in the 
SFR-age-metallicity-reddening space. Figure 5 displays a 
comparison between the SEDs of H3 and H5. The wave- 
length interval sampled by the broad-band filters in the 
restframe of H3 and H5 goes from the Lyman 912 A break 
to ~ 4400 A. The region at 1500 A is sampled by the 
Iw band, and the bands from R to J give the slope of 
the UV continuum, whereas the K' band is sensitive to 
the Balmer and 4000 A breaks (roughly the restframe B). 
Thus, we have in principle a good database to roughly 
constrain the stellar population in these objects. When 
the IMF and the upper mass limit for star-formation are 
fixed, the allowed parameter space can be roughly con- 
strained. The presence of Lya in emission points towards 
a star-forming system. Again, we used the GISSEL98 code 
for this exercise, taking into account that the two galax- 
ies are necessarily dominated by massive OB stars at the 



Restframe Wavelength (A) 




Restframe Wavelength (A) 




Fig. 5. Spectral energy distribution of H3 (top) and H5 (bot- 
tom), as determined from broad-band photometry. Fluxes are 
arbitrarily normalized to the mean flux in the Iw filter, and we 
use the same flux-scale for both sources. The best fit constant 
star-formation and burst models are displayed by thick and 
thin lines respectively. For H3, the best constant SFR has 0.02 
Zq metallicity, A v = 0.35 and age 2.4 Gyr (\ 2 = 0.85), and a 
burst model has 0.2 Zq metallicity, A v = 0.60 and age 0.004 
Gyr (x 2 = 0.91). For H5, the equivalent values are 2.5Zq, age 
0.51 Gyr (\ 2 = 0.83) for the constant SFR, and Zq, age 0.012 
Gyr (x 2 = 0.99) for the burst, all of them with A v = 0.0. See 
text for more details. 
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Figures 6a and 6b show the likelihood maps corre- 
sponding to H3 and H5 respectively, at z=4.05. For each 
SFR model, we have computed a reduced x 2 for 220 dif- 
ferent ages of the stellar population, 61 extinction values 
(ranging between A v = 0.0 and 3.0 magnitudes), and 5 dif- 
ferent metallicities (0.02Z Q , 0.2Z Q , Z & = 0.02, 2.5Z Q and 
5Zq). The most probable regions in this parameter-space 
are displayed in dark (permitted regions). The scaling di- 
rectly corresponds to the confidence level as derived from 
the x 2 value. The shaded regions enclose the 3a contours 
(confidence level of 99%). The likelihood maps were com- 
puted using two independent sets of filters: g,V w ,H,I w ,3,K 
and g,V w ,r,I w ,3,K, and the averaged intergalactic absorp- 
tion in the Lyman forest (see Sect. 3.2). In all the cases 
we use I w instead of I because the seeing and sampling 
are better in this filter. The results obtained with each 
independent set of filters are qualitatively the same, and 
the two la regions in the parameter space coincide. Each 
point on the likelihood maps presented in Figures 6a and 
6b corresponds to the most restrictive value obtained from 
the two independent sets (i.e. the lowest likelihood value). 
We also discuss in each particular case the results ob- 
tained when using only the spectral region redwards of 
Lya (hereafter rcdSED), a region which is not affected by 
uncertainties on the intergalactic absorption. In order to 
retain good resolution at short time scales, a logarithmic 
scale is used to display the age of the stellar population. 

According to Figures 6a and 6b, the stellar population 
in H5 is better constrained than in H3 in terms of A v . 
The reason for this is probably the relatively redder spec- 
trum of H3 compared to H5. The weight of the old stellar 
population is more important in H3 than in H5 and, as a 
consequence, we expect an enlarged set of permitted re- 
gions on the reddening-age plane. In a short-burst model 
without reddening, H3 should be older than H5. Figures 
6a and 6b also display the degeneracy in metallicity ver- 
sus reddening for the two objects. H3 and H5 do not seem 
to be highly reddened: the maximum restframe A v ever 
attained at a 3a level is ~ 1.1 magnitudes for H3, and it 
remains below ~ 0.8 magnitudes at la. This is not surpris- 
ing given the redshift of these objects and the photometric 
selection. The age of the stellar population in the constant 
star-formation models is unconstrained at 3a level in all 
the cases. Table 3 summarizes the permitted domains in 
the parameter space for H3 and H5. 

H3 is well fitted by both the burst and the con- 
stant star-formation models. Within the 99% confidence 
level, the 5Zq metallicity is excluded, whatever the star- 
formation model, the reddening or the age of the stellar 
population. All the solutions are compatible at la with 
a moderate A v ~ 0.4 to 0.6 magnitudes. The best fit at 
la with a burst model is obtained for metallicities below 
solar: 0.2 Z Q metallicity, A v = 0.60 and age 0.004 Gyr 
(X 2 = 0.91) At 3cr, the maximum age for the burst model 
is 0.1 Gyr, and the maximum value for A v is 1.1 magni- 
tudes. When we use only the redSED, the 3a domain is 



enlarged, and the maximum age and A v permitted for a 
burst model are 0.29 Gyr and 1.6 magnitudes respectively. 
The best fits of H3 are given by a constant star-formation 
model, 0.02 Zq metallicity, with A v = 0.35 and age 2.4 
Gyr (x 2 — 0.85), but this age is in fact incompatible with 
the age of the universe (1.17(2.16) Gyr with the adopted 
cosmology and go=0-5(0.1)). This is not a problem because 
the la domain for this model encloses solutions up to Zq 
metallicity, with ages starting at 0.18 Gyr and limited by 
the age of the universe. A v ranges between 0.1 and 1.1 
magnitudes at 3a, and it increases up to 1.6 magnitudes 
when using the redSED. 

H5 is well fitted by the burst and the constant star- 
formation models, provided that the A v value keeps below 
0.15 magnitudes typically. The best fit of H5 is given by 
a constant star-formation model of 2.5 Zq, age 0.51 Gyr 
and A v = 0.0 (x 2 = 0.83). The best fit by the burst model 
has Z Q , age 0.012 Gyr and A v = 0.0 (x 2 = 0.99). For 
the burst model, the age and the reddening are well con- 
strained. All the solutions within 2a are compatible with 
very small values of A v (A v < 0.3), whatever the model or 
the age of the stellar population. This result is also found 
when fitting only the redSED. At 3a, the maximum age 
of the burst is 0.07 Gyr, and up to 0.14 Gyr whith the 
redSED. The best fit models have solar metallicities or 
higher than solar. For the constant star-formation model, 
the reddening is well constrained. The ages are also well 
constrained (between 0.07 and the age of the universe), 
provided that we can exclude the highest metallicity. 

The best-fit SEDs were used to compute the SFR val- 
ues and the absolute magnitudes involved in the bursts 
from the observed magnitudes and mean fluxes. The lens- 
corrected luminosity observed within the I\y band is 
0.6(1. 7)10 45 /i 50 2 erg s" 1 for H3, and 2.6(7.4) lO 44 ^ 2 erg 
s _1 for H5, with go = 0.5(0.1), assuming A v = 0. In the 
case of H3, the best fit model has A v ~ 0.5, thus the cor- 
rected luminosity in this filter (restframe ~ 1450 - 1700 A) 
is in fact 2.8 times higher. When using the best-fit mod- 
els to scale the fluxes (Fig. 5), we obtain L(1500 A)= 1.24 
(3.5) lO 42 /^ 2 erg s" 1 A" 1 for H3 (corrected for A v = 0.5), 
with q = 0.5(0.1), and L(1500 A)= 1.9 (5.4) lO 41 ^ 2 erg 
s _1 A -1 for H5 (A v = 0). These values are not strongly 
dependent on the model SED used for the scaling. The 
SFRs derived from the best-fit star-forming models, scaled 
to the observed L(1500 A), are 311(110)/i^ 2 M Q /yr for 
H3 (corrected for A v = 0.5) and 54(19)/i 50 2 M Q /yr for H5 
(A v = 0), and only ~ 1/3 of these values correspond to 
stars with m> 1M . The main uncertainty is the A v for 
H3, which could set the SFR to values ~ 3 times lower 
or higher than the best solution quoted above. The ab- 
solute B magnitudes were derived from the observed K 
magnitudes (restframe ~ 3900 - 4430 A), using the best-fit 
models for a detailed scaling: Mb — —22.9 (—24.1) and 
M B = -20.7 (-21.9) for H3 and H5 respectively, with 
<7o = 0.5(0.1). In the case of H5, this result is almost inde- 
pendent on the model SED because A v is always small; in 
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Fig. 6. Likelihood map of H3 (left) and H5 (right) showing in dark the most probable regions and the degeneracy in the 
parameter space defined by SFR type, age, metallicity and reddening. Dotted and solid lines give the age limit corresponding 
to z=4.05, with go = 0.5 and 0.1 respectively. Metallicities are given in units of log(Z)=log(Z/i?0), and the solar metallicity 
is Zq = 0.02. The shaded regions enclose the 3<r contours (confidence level of 99%, or a likelihood value of 1%, according the 
scaling displayed at the top). A logarithmic scale is used for ages in order to retain good resolution at short time scales. 



Table 3. Permitted domains in the parameter space of H3 and H5. The values given are rough limits, the correlation between 
the different parameters is displayed in figures 6 and 7. 





Model 


lcr 








3cr 










Age (Gyr) 


A v (mags) 


Metallicity (Z e ) 




Age (Gyr) 


A v (mags) 


Metallicity (Z @ ) 


H3 


Burst 


0.003 - 0.007 


0.5 - 0.8 


0.02 - 0.2 




0.0 - 0.10 


0.0 - 1.1 


0.02 - 2.5 


H3 


Constant SFR 


0.18 - 1.17(2.16) 


0.2 - 0.6 


0.02 - 1.0 


0.0 


- 1.17(2.16) 


0.1 - 1.1 


0.02 - 2.5 


H5 


Burst 


0.003 - 0.016 


0.0 - 0.25 


0.2 - 5.0 




0.0 - 0.07 


0.0 - 0.7 


all 


H5 


Constant SFR 


0.0001- 1.17(2.16) 


0.0 - 0.15 


all 


0.0 


- 1.17(2.16) 


0.0 - 0.6 


all 



the case of H3, the 3cr uncertainties in A v translate into 
~ ±0.7 magnitudes of uncertainty in Mg. 

The shape of the continuum at wavelengths shorter 
than ~ 8000 A (restframe ~ 1500 A), up to the I band, 
is relatively insensitive to the metallicity for metallicities 
higher than solar. The likelihood map is almost insensitive 
to age for ages below ~ 10 6 years, as expected given the 
sampling in stellar masses in the evolutionary tracks used 
by GISSEL98. Nevertheless, such time-scales seem irrele- 
vant here, the best solutions being older than this limit. 
These results also are generally insensitive to the choice of 
the IMF. Nevertheless, the permitted regions in the like- 



lihood maps show some dependency on the upper mass 
limit assumed for star-formation. In particular, when this 
limit is set to a value as low as 10M Q , the permitted region 
for burst models is shifted towards younger age values for 
the stellar population, irrespective of the metallicity. This 
change on model details has very small influence on the 
SFRs derived above. 

5. Discussion and conclusions 

The two sources H3 and H5 are at the same redshift, but 
their stellar populations are noticeably different. Both of 
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them are well fitted by burst or continuous star-formation 
models. Metallicities much higher than 2.5Zq seem to be 
excluded by the photometric data in the case of H3, and 
they have a lower probability in H5. Both galaxies are 
roughly compatible with solar metallicities within the la 
confidence level: Z < Z Q for H3 and Z > 0.2Z & for H5. The 
stellar population seen in H5 is younger than in H3, on the 
basis of a burst model with A v = 0. The main source of un- 
certainty when deriving the restframe properties of these 
objects is the value of A v , and this is clearly shown in the 
case of H3 where the 3cr uncertainties are ~ ±0.5 mag- 
nitudes. The two sources are intrinsically bright: roughly 
M B for H5 and M* B - 2 (±0.7) magnitudes for H3 (here 
M B is the local value, from Loveday et al. 1992). This, 
combined with the high gravitational amplification and 
the presence of relatively strong emission lines, has allowed 
us to obtain a spectroscopic redshift for these sources using 
a 4m telescope. Among the high redshift candidates in our 
sample behind A2390, 4 additional ones are at the same 
photometric redshift. A subsequent spectroscopic survey 
using a 8m class telescope is needed to go further on this 
study, especially to estimate more precisely the metallici- 
ties from UV absorption lines. 

When comparing our results on H3 to those by Frye & 
Broadhurst (1998) and Bunker et al. (1998), we find them 
in fairly good agreement, after correction for the difference 
in the amplification factor, which is 3 times higher in their 
lens model. The age of the system and the reddening value 
given by Bunker et al. for H3 are included within our 2a 
best solution region, which is a remarkable result taking 
into account that both the photometric data (HST ex- 
cepted) and the analysis are completely independent. The 
high amplification factor obtained by Frye & Broadhurst 
(1998) corresponds to the maximum value attained in our 
model for H3, on the region neighbouring the critical line, 
but the surface-averaged value in our case is 3 times lower. 

Compared to other known z > 3 galaxies, H3 and H5 
belong to the bright end of the field population (Lowen- 
thal et al. 1997; Steidel et al. 1996a and 1996b). They are 
slightly brighter than the z ~ 4 objects found by Trager et 
al. (1997) in the cluster lens C10939+4713, although the 
amplification factor is poorly known in this case. They 
are also intrinsically brighter than the z = 5.34 galaxy 
found by Dey et al. (1998). These sources are not resolved 
in their width, where the lens inversion is limited by the 
resolution of the WFPC2 images (about 0.1", 0.9(1. 6)h^ 
kpc with (jo = 0.5(0.1)), but they are resolved on their 
length. H3 and H5 are splitted into several small, com- 
pact and bright subclumps, all aligned towards the same 
direction. The total length of these emitting regions is sim- 
ilar to that of the compact cores of the 2.5<z<3.5 field 
sample by Steidel et al. (1996a). H3 is more elongated 
than H5, and probably more dusty. The linear separation 
between H3 and H5 (150(252)^^ with q = 0.5(0.1)) 
and their peculiar morphologies strongly point towards a 
hierarchical merging process as a likely scenario for the 



formation of the brightest spheroids. Their photometric 
SEDs are compatible with a wide range of metallicities at 
a 3a level, all with Z > 0.02%Z Q . Nevertheless, as shown in 
Section 4, broad band photometry does not allow a precise 
estimate of the metallicity. A spectroscopic study using 
UV restframe absorption lines, or near-IR spectroscopy, 
is urgently needed for these purposes. A relatively high 
metallicity for the bright subclumps would imply that we 
are actually seeing an advanced step in this merging pro- 
cess, where stars are forming from a metal enriched gas 
(see also Trager et al. 1997; Lowenthal et al. 1997, Baugh 
et al 1997; Moscardini et al. 1997). It is worth noting that 
we are dealing with local conditions in these star form- 
ing systems. The light detected is mainly emitted in very 
small and compact regions (~ l/i^~ kpc), which could be 
highly enriched compared to the remainder of the source. 
This is fully compatible with the best fit metallicities being 
> 0.02%Z Q in these systems, because in this case the star 
formation becomes a local and efficient process, where the 
cooling rate is enhanced by metallic atoms allowing the 
formation of molecules and increasing the dust opacity. 
Such a process is discussed in details in a recent paper by 
Spaans & Carollo (1998). 

The uncertainty in the amplification factor for H3 and 
H5 is 0.3 magnitudes. This means that the intrinsic lumi- 
nosities and SFRs are known with an accuracy of ~ 30%. 
This source of error has the same importance than the 
model uncertainties for a relatively well constrained SED 
(i.e., the uncertainty of ~ 0.5 magnitudes in A v for H3). 
It is worth noting that H3 and H5 are multiple images, 
well modelled compared to other images in this or other 
cluster lenses. This gives an idea of the limitation aris- 
ing from lens modelling when using clusters as gravita- 
tional telescopes to access the background sources. Only 
the well constrained clusters are actually useful for this 
programme. 

H3 and H5 are the first spectroscopically confirmed im- 
ages of sources at z > 2 in this cluster, a redshift domain 
which is well constrained by the set of filters used here. 
The selection of high-redshift candidates using a photo- 
metric redshift approach, including the near-IR bands, is 
strongly supported by the present results. For most statis- 
tical purposes, photometric redshifts should be accurate 
enough to discuss the properties of these extremely dis- 
tant galaxies. Conversely, the spectroscopic confirmation 
of the redshifts of such gravitationally amplified sources 
could help on the calibration and improvement of the pho- 
tometric redshifts techniques up to much fainter limits in 
magnitude compared to field surveys. 
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